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ABSTRACT 

Aims. We present the most precise light curve ever obtained of a detached eclipsing binary star and use it investigate the inclusion of 
non-linear limb darkening laws in light curve models of eclipsing binaries. This light curve, of the bright eclipsing system f3 Aurigae, 
was obtained using the star tracker aboard the wire satellite and contains 30 000 datapoints with a point-to-point scatter of 0.3 mmag. 
Methods. We analyse the wire light curve using a version of the ebop code modified to include non-linear limb darkening laws and 
to directly incorporate observed times of minimum light and spectroscopic light ratios into the photometric solution as individual 
observations. We also analyse the dataset with the Wilson-Devinney code to ensure that the two models give consistent results. 
Results, ebop is able to provide an excellent fit to the high-precision wire data. Whilst the fractional radii of the stars are only defined 
to a precision of 5% by this light curve, including an accurate published spectroscopic light ratio improves this dramatically to 0.5%. 
Using non-linear limb darkening improves the quality of the fit significantly compared to the linear law and causes the measured radii 
to increase by 0.4%. It is possible to derive all of the limb darkening coefficients from the light curve, although they are strongly 
correlated with each other. The fitted coefficients agree with theoretical predictions to within their fairly large error estimates. We 
were able to obtain a reasonably good fit to the data using the Wilson-Devinney code, but only using the highest available integration 
accuracy and by iterating for a long time. Bolometric albedos of 0.6 were found, which are appropriate to convective rather than 
radiative envelopes. 

Conclusions. The radii and masses of the components of ySAur are R A = 2.762 ± 0.017 R Q , R B = 2.568 ± 0.017 R Q , M A = 2.376 ± 
0.027 M Q and M B = 2.291 ± 0.027 M Q , where A and B denote the primary and secondary star, respectively. Theoretical stellar 
evolutionary models can match these parameters for a solar metal abundance and an age of 450-500 Myr. The Hipparcos trigonometric 
parallax and an interferometrically-derived orbital parallax give distances to f3 Aur which are in excellent agreement with each other 
and with distances derived using surface brightness relations and several sets of empirical and theoretical bolometric corrections. 

Key words, stars: fundamental parameters — stars: binaries: eclipsing — stars: binaries: spectroscopic — stars: distances — stars: 
evolution — stars: individual: beta Aurigae 



5h I 1. Introduction but the coolant for the main camera was lost only a few days 

after launch. Since then the star tracker on wire has been used 

The study of detached eclipsing binaries (dEBs) is our primary t0 monitor the variability of brigh t stars of all spectral types, as 

source of accurate physical properties of normal stars ( |Andersen| desC ribed in the review by Bruntt ( 2007). The main aims of the 

\1993- By modelling their light curves and deriving radial veloc- current work are t0 accurate i y measure the physical properties 

ities of the two components from their spectra, it is possible to of the systems wd t0 investigate whether the currently available 

measure the masses and radii of the two components to accura- eclipsing binary light curve models can provide a good fit to 

cies of better than 1 % (e.g.,| Southworth et al. || 2005ct | Lacy et al. | observational data with measurement uncertainties substantially 

|2005|). The effective temperatures and luminosities of the stars below j mmag We expect the answer to the second quest ion to 

can be found by a variety of different methods, including pho- be < yes < ; but this matter is bec0 ming increasingly important in 

tometric calibration, spectral analysis and measuring trigone- light of current an d forthcomin g space missio ns, such as Kepler 

metric parallaxes. These data can be used to measure accurate dB as ri et alJl2005h and CoRoT dMichel et alJl2006n . which aim 

distances to stellar systems ( Ribas et al. || 2005 |; | Southworth et al. | t0 obtain extremely high-quality light curves of a large number 

2005a) and to critically investigate the accuracy of the predic - Q f stars 
tions resulting from theoretical stellar models (e.g., Ribas 2006). 

We are undertaking a project to obtain accurate and extensive Our first paper in this series presented the serendipitous 

light curves of selected bright dEBs using the Wide Field Infra- discovery that the bright star i/f Centauri is a long-period dEB 

red Explorer (wire) satellite, wire was launched in March 1999 dBruntt et all 12006, hereafter Paper I). The wire light curve 

of this object shows deep total eclipses arising from a sys- 

Send offprint requests to: John Southworth tern containing two stars of quite different sizes. dEBs con- 

* Table 2 is only available in electronic form at the CDS via taining two dissimilar stars are particularly useful because 

anonymous ftp to cdsa rc . u-strasbg . fr ( 130.79.128.5) or via the co-evolutionary nature of the two components can put 

http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ strong constraints on the predictions of theoretical models 



2 



Southworth, Bruntt & Buzasi: WIRE satellite photometry of /? Aurigae 



(Southwo rth et al.ll2004bl) . However, one of the most interesting 
aspects of this work is that we were able to measure the relative 
radii of the two stars (radii expressed as a fraction of the orbital 
size) to a precision of better than 0.2%, where the uncertainties 
were estimated from extensive Monte Carlo simulations but did 
not include error contributions from some sources such as the 
use of different algebraic laws to describe the limb darkening of 
the stars. 

pAurigae 

In this work we present a wire light curve for the dEB /? Aurigae 
(Menkalinan), which is composed of two Al IV stars orbiting 
each other with a period of 3.96 days (Table[TJ. This light curve 
contains 30015 datapoints with a point to point scatter of less 
than 0.3mmag. We chose /3Aur as a target for the wire satel- 
lite due to its brightness, which not only made our observations 
more precise but also makes it a very difficult object to study 
with modern telescopes. In addition, it has an accurate trigono- 
metric parallax from the Hipparcos satellite which means the 
effective temperatures and surface brightnesses of the stars can 
be measured directly and used to calibrate these quantities for 
other objects. 

P Aur was one of th e first stars dis covered to be a spectro- 
scopic bin ary system, by Maury ( 1890), and an orbit was calcu- 
lated by iRambaua (11891) from the published observations. It was 
subsequently discovered bv lStebbinsI d 1 9 1 lb to undergo shallow 
eclipses, and so was also one of the first known eclipsing bina- 
ries. A detailed discussion of early work, and an investigation 
into whet her the speed of light is dependent on wavelength, was 
given bv lBakerl(ll910l) . 

The brightness of /3 Aur and the shallowness of its eclipses 
has made it difficult t o ob serve photometrically, and only 
the study by iJohanserJ dl97ll) has provided accurate results. 
Johansen obtained light curves of ySAur in seven passbands, 
four of which closely resemble t he Stromgren uvby system , and 
solved them using the method of Russell & MerriJ d 1952b . 

However, the brightness of ySAur and the relatively low 
rotational velocities of the two components means that it is 
observationally one of the most straightforward spectroscopic 
binaries to st udy. Spec t roscop i c radia l veloc i ties ha v e bee n 
presented by iRambautl d 189 lb . IVogell (119041). ISmithl d!948b . 
iGaleotti & Guerrero! d 19681) and [Popper & Carlosl d 19701) . the 
most accurate and detailed being that by Smith. In addition 
to thes e work s, the spectrum of /3Aur h as be en consider ed by 
iMaurvl d 1898b . IStruve & Driscoill d 1953b a nd iTovl d 19691) . the 
most recent study (Lvubimk ov et al.| [l 996) demonstrating that 
the two components have almost identical spectra and exhibit 
en hanced metallic l ines t ypical of the Am class of stars. 

iHummel et alj d 1995b have interferometrically resolved the 
orbit of ySAur using the Mark III Optical Interferometer. This 
means that the distance to the system can be obtained straightfor- 
wardly by comparing their angular orbital semimajor axis to the 
linear value measured using spectroscopy. In addition, Hummel 
et al. have measured the orbital inclination, which can be com- 
pared with the value derived from the light curve to ensure con- 
sist ency. 

Nordstrom & Johansen ( 1994, hereafter NJ94) derived accu- 
ra te absolute dime nsions of /? Aur by analysing the light curves 
of |Johansenldl97lb with the eb op code and adopting the spectro- 
scopic orbits from lSmiThl d!948h . It is well known that for dEBs 
containing two stars of almost the same size, particularly if they 
exhibit shallow eclipses, the ratio of the ste llar radii is poorly de- 
termined by the light curves ('Pop per! 1984b . so NJ94 constrained 



Table 1. Identifications and astrophysical da ta fo r fi Au r. 
References: (1) iHoffleit & Jaschekl dl991b: (2) 
ICannon & Pickering! (fl918b; (3) iPerrvman et al I dl997b : (4) 



NJ94: r5^Lvubimko vetaiTdT9^ :76T lH0g et al.ldl997b " 





/3Aur 


References 


Bright Star Catalogue 


HR 2088 


1 


Henry Draper Catalogue 


HD 40183 


2 


Hipparcos Catalogue 


HIP 28360 


3 


^2000 


05 h 59 m 31.7229 s 


3 




+44° 56' 50.758" 


3 


Spectral type 


AlmIV + Aim IV 


4.5 


Trigonometric parallax 


39.72 ± 0.78 mas 


3 


Distance 


25.18 ± 0.49 pc 


3 


Tycho B magnitude (B T ) 


1.989 ± 0.005 


6 


Tycho V magnitude (W) 


1.904 ± 0.005 


6 



their light curve solution with an accurate spectroscopic light ra- 
tio. The resulting properties of the stars have s ince been used for 
studies of the helium-metal enrichment law (Ribas et al. 2000) 
and th e fundamental effective temperature scale dSmalley et alJ 
12002b . 



2. Observations and data reduction 

/3 Aur was observed with the wire star tracker from 2006 March 
30 to April 20. During each 93.2 min orbit of Earth wire 
switched between two targets in order to minimise the influ- 
ence of scattered light from the daytime side of Earth. In the 
data reduction we discarded the datapoints at the very begin- 
ning and the end of each orbit, resulting in a duty cycle of 25%. 
Observations were obtained during e very orbit for 21 da ys. The 
wire passband is approximately V+R dBruntt et al.l2007b . It may 
in future be possible to obtain a more precise description of the 
response function by considering the the total set of observations 
obtained using the wire star tracker. 

The raw data comprises about 950000 CCD windows of 
8x8 pixels collected at a rate of 2 Hz. To extract th e light curve 
we us ed the photometric pipeline described by Brunt retaTI 
d2005b . The background level was very low in these observa- 
tions (20 ADU or 300 e~), and hence the photometric precision 
is mainly determined from Poisson statistics. We binned the light 
curve to four datapoints per minute, giving a total of 30015 dat- 
apoints after rejecting a few clear outliers. 

The point-to-point scatter in this light curve is 0.28 mmag, of 
which the Poisson noise accounts for 0.21 mmag as each orig- 
inal observation contains about 770000 e~. We calculated the 
noise in the Fourier amplitude spectrum in the frequency range 
2-3 mHz. The mean amplitude level is (A) = 2.7 //mag, which 
corresponds to y/{A) * N/n = 0.27 mmag per data point. The 
obser ved noise leve l for P Aur is higher by 29% than the Poisson 
noise. iBruntt et~aT] (120051) analysed wire data of Procyon and 
their measured noise level was 15-20% higher than the Poisson 
noise. These higher noise levels may be due to instrumental ef- 
fects like satellite jitter combined with variable sub-pixel effi- 
ciency in the CCD. 

After subtracting a preliminary light curve fit we clearly saw 
long-period drifts at low amplitude. In particular, for the first 
day of observations the star appears to slowly decrease in bright- 
ness by 2 mmag and during the last two days of observations the 
star apparently again became brighter by 2 mmag. This change is 
tightly correlated with a slight increase in the background level, 
and is clearly not intrinsic to /3 Aur as it is not affected by the 



Southworth, Bruntt & Buzasi: WIRE satellite photometry of /3 Aurigae 



3 



1 .90 



„ ,, ,i „ 

, '" , / ' , 



1 .92 



1 .94 



£ 1.96 I 



1 .98 



53825 53830 53835 53840 53845 

HJD - 2400000 

Fig. 1. The wire light curve of ft Aur. Each group of points represents continuous observation for a part of one orbit of the wire 
satellite around Earth. The overall magnitude has been adjusted to match the Hipparcos H p light curve but is not exact because the 
passband of wire is unknown. 



presence of eclipses. In the Fourier spectrum of the residual light 
curve we detected two peaks at f\ = 1.020 and = 5.148 
cycle d~', both with amplitudes of just 70//mag. The former is 
likely related to a slight change of the phase of the illuminated 
Earth. The second frequency corresponds to exactly one third of 
the orbital frequency of WIRE (/wire = 1 5. 449 cycle d~'). We 
cannot be sure why this is visible in the data. To remove the ef- 
fects mentioned above from the light curve, we used a high-pass 
filter on the observations. This was done on the residuals of the 
preliminary light curve fit, and four iterations were made to en- 
sure the depths or shapes of the eclipses were not affected. The 
high-pass filter has a cut-off which removes drifts in the data 
corresponding to periods longer than 24 hours. 

We note that in the analysis of WIRE data of the binary 
i^Cen (Paper I) we found two significant peaks at f\ = 1.996(2) 
and fi = 5.127(3) cycle d , both with amplitudes around 
0.2mmag. In Paper I we proposed that these frequencies were 
due to pulsations intrinsic to the star. However, at the time of 
those observations the orbital frequency of WIRE was slightly 
lower at /wire = 15.348 + 0.001 cycle d , which is exactly one 
third of fi as we also find here for ft Aur. It is therefore very 
likely that the variation found at the frequency |/wire is not due 
to variations intrinsic to either ift Cen or ft Aur. 



3. Light curve analysis 

The wire light curve of ft Aur is shown in Fig.Q] It lasts for just 
under 21 days and contains 11 eclipses. The observations are 
presented in Table|2] (available in its entirety in the electronic 
version of this work). We define phase zero to occur at the mid- 
point of the primary eclipse, which is deeper than the secondary 
eclipse. We will refer to the eclipsed component at this point as 
star A and the (cooler) eclipsing body as star B. For ft Aur star A 
is slightly larger and more massive than star B, as well as being a 
little hotter. The phased wire light curve is shown in Fig.|2]with 
the residuals of the best fitting model (see below). 



Table 2. wire differential-magnitude light curve for ft Aur. The 
complete table is available in the electronic version of this work. 



HJD 


Differential magnitude 


2453825.314879 


0.0195836 


2453825.315059 


0.0197159 


2453825.315238 


0.0193520 


2453825.315417 


0.0190368 


2453825.315597 


0.0192776 


2453825.315776 


0.0186899 


2453825.315956 


0.0187692 


2453825.316135 


0.0190573 


2453825.316314 


0.0187044 


2453825.316494 


0.0184744 



For th e main analysis of the w ire light curve we have used 
jkteboi{3 ( Sout hworth et alJl2004bj). wh ich is based on the ebop 
code dPopper & Etzelll 198 ltlEtzellf 1 98 lb , which in turn uses the 
NDE biaxial ellipsoidal model for well-d etached eclipsing bi- 
naries dNelson & Davislll972tlEtzeiT975T) . jktebop incorporates 
several modifications to the standard code, particularly in in- 
put/output and the use of t he Levenberg-M arquardt minimisa- 
tion algorithm (mrqmin; see lPress et al.ll 19921 p. 678) for finding 
the optimal model parameters for describing a dataset. The rel- 
ative radii of the stars are defined to be r& — ^ and — ^f, 
where and Rr are the absolute radii of the stars and a is the 
orbital semimajo r axis, jktebop solves for the sum and ratio of 
the relative radii (South worth et al.ll2005d) as these are more di- 
rectly related to the morphology of dEB light curves than the 
individual radii. 

For this study we have implemented three additional modifi- 
cations: the incorporation of non-linear limb darkening and the 
direct inclusion of spectroscopic light ratios and times of mini- 
mum light as observational data. These will be discussed below. 



1 jktebop is written in fortran77 and the source code is available at 
http : //www. astro . keele . ac . uk/~ jkt/ 
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Fig. 2. Phased wire light curve of /3 Aur. The residuals of the best fit are shown in the lower panel. The best fit has not been plotted 
in the upper panel as it completely overlaps the data (as can be seen from the residuals). 



3.1. Period analysis 

A common method of determining the linear orbital ephemeris 
of an eclipsing binary is to collect times of minimum light, 
spanning as large a time interval as possi ble, and fit them 
with a straight line against cycle number (e.g JSouthworth et al.l 
2005c). This procedure is not a good approach here due to the 
properties of the wire light curve: it is possible to obtain individ- 
ual times of minimum by fitting Gaussian functions separately 
to each eclipse, but these contain systematic errors in the mea- 
sured times as the eclipse shape is slightly non-Gaussian and 
individual eclipses are poorly sampled. An improvement over 
this would be to fit each eclipse with a light curve model, but the 
best way to measure an orbital period for an eclipsing binary is 
to fit for it directly using the whole light curve. This is often not 
possible because the available times of minimum light are not 
accompanied by the data from which they were measured. 

We have therefore modified jktebop to treat times of min- 
imum light themselves as observational data and include these 
alongside the light curve when finding the best fit. The inclusion 
of disparate types of data obviously requires the measurement 
errors of each type to be acc urately determ ined so their relative 
weightings are correct (see IWilsonlll993l) . We have measured 
times of minimum light from the wire light curve using Gaussian 
functions. We report these for the interested reader (Tabled but 
have not used them in our analysis. 

A detailed perio d analysis for /3Am was given by 

lJohansen~& S0rensen (1997) on the basis of several observed 
times of minimum light and spectroscopically or interferomet- 
rically determined times of conjunction. Whilst there have been 
occa sional sugg estions that the period of /3 Aur is not constant 
(see IVogel|[l904l) . Johansen & S0rensen found no evidence for 
a variable orbital period. We have chosen to include only those 



times of minimum light which were obtained via photoelectric 
or CCD photometry as these measurements are more accurate 
than the times of conjunction obtained in other ways. 

In order to derive the orbital period of /3 Aur without acciden- 
tally introducing cycle count errors we began with a preliminary 
light curve fit to the the wire observations. Once we had deter- 
mined the best period for this datas et we added in the Hipparcos 
light curve dPerrvman et alj 119971) . with appropriate observa- 
tional errors, to give a more accurate period (see also Paper I). 
The slightly different bandpasses for these two datasets is not 
important because the quality and quantity of the wire obser- 
vations are both orders of magnitude better than the Hipparcos 
data. Inclusion of the latter makes a negligible difference to the 
parameters of the fit, except for increasing the accuracy of the 
orbital period determination. 

We then introduced th e ti mes of minim um light from 
iJohansen & S0rensenl (1 19971) and I Johansen! dl97ll) before finally 
extending the dataset to cover almost a century with the time of 
minimum from Steb binsl d 1 9 1 ll) quoted by Johansen & S0rensen. 
The quality of the fit to the times of minimum is excellent 
(Table|3j, and from consideration of these data and the other 
times of minimum used by Johansen & S0rensen we can find 
no evidence for a change in period. 



3.2. Direct inclusion of a spectroscopic light ratio in the light 
curve fitting process 

For dEBs which contain two similar stars undergoing partial 
eclipses, the sum of the radii (r& + r%) is well defined as it 
depends mainly on the orbital inclination and duration of the 
eclipses. The ratio of the radii (k = ^) can be very poorly 
defined by the observed light curves of such systems because 
quite different values can give very similar eclipse shapes (e.g. 
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Table 3. Times of minimum light for /3 Aur compared to the calculated values for our chosen model fit. Times are given in HJD - 
2400000. 

* Note that the times measured from the wire data are expected to be affected by systematic errors (see Section lXTl ) and were not 
included in the fit. They are repeated here for completeness and for the interested reader. 



Orbital cycle 


Observed time 


Calculated time 


O-C 


Reference 


-8790.0 


19018.3848 


± 0.0046 


19018.38497 


-0.00017 


Stebbins (1911) 


-4289.0 


36842.5537 


±0.0013 


36842.55528 


-0.00158 


Johansen (1971) 


-4085.0 


37650.4025 


±0.0011 


37650.40481 


-0.00231 


Johansen (1971) 


-3719.0 


39099.785 


±0.0011 


39099.78191 


+0.00309 


Johansen (1971) 


-1030.0 


49748.353 


±0.0023 


49748.34755 


+0.00545 


Johansen & S0rensen (1997) 


-855.0 


50441.3552 


±0.0035 


50441.35573 


-0.00053 


Johansen & S0rensen (1997) 


0.0* 


53827.19586 


±0.00045 


53827.19567* 


+0.00019 


This work 


0.5* 


53829.17566 


± 0.00047 


53829.17570* 


-0.00004 


This work 


1.0* 


53831.15541 


± 0.00054 


53831.15572* 


-0.00031 


This work 


1.5* 


53833.13593 


±0.00029 


53833.13574* 


+0.00019 


This work 


2.0* 


53835.11562 


±0.00032 


53835.11577* 


-0.00015 


This work 


2.5* 


53837.09603 


± 0.00042 


53837.09579* 


+0.00024 


This work 


3.0* 


53839.07576 


±0.00019 


53839.07581 * 


-0.00005 


This work 


3.5* 


53841.05571 


± 0.00047 


53841.05584* 


-0.00013 


This work 


4.0* 


53843.03598 


± 0.00024 


53843.03586* 


+0.00012 


This work 


4.5* 


53845.01563 


±0.00037 


53845.01588* 


-0.00025 


This work 



lAndersen et al.fl980l) . ft Aur is unfortunately an excellent exam- 
ple of this situation, despite the quantity and quality of the wire 
light curve. However, the light ratio of dEBs is normally very 
strongly dependent on k so an observed spectroscopic light ra- 
tio can be use d to constrain k and therefore the individual stel- 
lar ra dii (e.g. ISouthworth et all l2004at ISouthworth & Clausen! 
120071) . 

The usual way this constraint is applied is to map out how the 
light ratio depends on k and use the independent sp ectroscopic 
light ratio to specify the allowable values of k (e.g. lTorres et alj 
2000). However, a more rigorous procedure is to include the 
spectroscopic light ratio directly as an observation to be used 
during the minimisation process. This also greatly simplifies the 
error analysis, jktebop has been modified to explicitly include 
this type of observable quantity in the fitting process. As with 
inclusion of the times of minimum light, it is important to en- 
sure that accurate measurement errors are adopted when mixing 
different types of data in one fitting process. 

In preliminary fits to the wire light curve we found that r& 
and re were determined to a precision of only 5% (uncertainty 
from Monte Carlo simulations; see below). This was worse than 
expected, and quite distant from the target of 2% w hich is the 
minimum for obtaining use ful results for most dEBs (iAndersenl 
1 1 99 ll 1 1 998t iGimenezI 1992b . To solve this problem in their anal- 
ysis of j3 Aur, NJ94 obtained a high-quality spectroscopic light 
ratio of 0.855 ± 0.016, which they used to constrain k and so ob- 
tain accurate values for and r%. This light ratio is consistent 
with and much more accurate th an the magnitude diff erences 
of Am = 0.2 + 0.2 mag given b vlHummel et alj d 19951) and of 
Am = 0.1 3 mag given bv lPetrid(ll939h . When we included the 
NJ94 light ratio in our fit we found that the precision of the radii 
improved dramatically from 5% to 0.5%. k is determined almost 
entirely by the light ratio in this solution. 

A possible problem in this approach is that NJ94's light ratio 
was obtained using the Mg II spectral line at 448 1 A but the wire 
passband is redder than this so the light ratio may be different. 
However, the two stars have very similar effective temperatures 
so the systematic error caused by this effect should be negligible. 
We note that it is not possible to prove this because the transmis- 
sion curve of the wire star tracker is not precisely known. In the 



future we will obtain a spectroscopic light ratio at a wavelength 
closer to the effective wavelength of the wire passband. 

3.2.1. Non-linear limb darkening 

A complication in fitting the light curve is the inclusion of limb 
darke ning (LD). The ebop code uses the linear LD law dRusselll 
Il912l) . which normally provides an entirely acceptable fit to the 
light curves of dEBs but is known to be a woefully inadequate 
repre sentation of the pr edictions of stellar model atmospheres 
(e.g. IVan Hamme| [T993). T here is evidence from stud ies using 
the Wilson-Devinney code dWilson & Devinnevlll971l) that in- 
cluding bi-parametric laws can cause a slight improvement in 
the fit and a systematic increase of the order of 0-2% in t he mea- 
sured fractional stellar radii (Paper I; lLacv et aTl I2005I) . which 
is certainly large enough to be important in the present case. 
This can be understood because non-linear LD laws generally 
result in a lower fl ux towards the stellar limb compared to the 
linear law (see e.g. lClaretll2000al:IClaret & Hauschildtll2003l) . so 
a larger area is required to produce the same amount of light at 
the limb. 

jktebop has therefore been modified to allow both the use of 
several non-linear LD lawfl and the ability to optimise both co- 
efficients for each star. The i mplemented laws are the quadratic 
dKopall [l950T) . s quare-root dDiaz-Cordoves & Gimenez! 1 19921) 
and logarithmic (Klinglesmith & Sobieskil ll970h . 

It is common practise when fitting the light curves of eclips- 
ing binaries to use LD coefficients obtained from theoretical 
model stellar atmospheres, partly because the light curve shape 
is not strongly dependent on the LD coefficients and a good set 
of data is needed t o be able to fit for them directly (Popper 1984; 
iNorth et all 1 1 9971) . However, this means that the final solution 
depends on theoretical predictions, so it is better to fit for the 
coefficients when possible. 

We have performed light curve solutions for /3 Aur using the 
linear, square-root and quadratic LD laws. In each case we in- 
cluded the LD coefficients for both stars in the set of fitted pa- 



2 Non-linear LD law s have previously been incorporated into ebop by 
iDfaz-CordovesI (1990). We are grateful to Dr. A. Gimenez for sending 



fir; 

us the relevant lines of code. See Gimenez & Dfaz-Cor doves ( 1993]) and 
Gimenez & Ouintana ( 1992). 
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Table 4. Best-fitting photometric parameters. These solution have been obtained by fitting the wire and Hipparcos light curve as 
well as the spectroscopic light ratio of NJ94 and several times of minimum light. Times are expressed in HJD - 2 400 000. 



Results using individual LD coefficients for each star 

Linear LD 



law 



Square-root LD law 



Quadratic LD law 



Sum of the relative radii 


r\ + r B 


0.30174 


±0.00007 


0.30297 


±0.00036 


0.30290 


±0.00047 


Ratio of the radii 


k 


0.9289 


±0.0083 


0.9311 


± 0.0083 


0.9300 


± 0.0084 


Surface brightness ratio 


.1 


0.9883 


±0.0028 


0.9812 


±0.0091 


0.9856 


±0.0126 


Orbital inclination (deg.) 


i 


76.862 


±0.004 


76.809 


±0.015 


76.803 


± 0.020 


Eccentricity term 


ecosw 


-0.000017 


±0.000006 


-0.000016 


± 0.000006 


-0.000015 


± 0.000006 


Eccentricity term 


esmco 


0.00191 


±0.00032 


0.00173 


± 0.00036 


0.00183 


± 0.00035 


Orbital period (days) 




3.96004673 


±0.00000015 


3.96004673 


±0.00000015 


3.96004673 


±0.00000015 


Time of primary minimum 




52827.195695 


±0.000011 


52827.195693 


±0.000011 


52827.195693 


±0.000011 


Linear LD coefficient 


"A 


0.382 


±0.006 


-0.056 


±0.147 


0.210 


±0.064 


Linear LD coefficient 


m b 


0.379 ±0.006 


-0.177 


±0.148 


0.187 


±0.060 


Non-linear LD coefficient 


Va 






0.703 


± 0.232 


0.248 


±0.096 


Non-linear LD coefficient 


Vb 






0.881 


±0.236 


0.277 


±0.094 


Relative radius of star A 


>A 


0.15643 


±0.00065 


0.15689 


± 0.00068 


0.15694 


± 0.00070 


Relative radius of star B 


H5 


0.14531 


±0.00069 


0.14608 


± 0.00072 


0.14595 


± 0.00074 


Light ratio 


^b/^a 


0.853 


±0.016 


0.855 


±0.016 


0.854 


±0.016 


Orbital eccentricity 


e 


0.00191 


±0.00032 


0.00173 


± 0.00036 


0.00183 


± 0.00035 


Periastron longitude (deg.) 


CO 


90.51 


±0.21 


90.51 


±0.25 


90.47 


±0.23 


Reduced x 1 


2 

X red 


1.332 


1.328 


1.327 



Results using the same LD coefficients for the two stars 

Linear LD law Square-root LD law Quadratic LD law 



Sum of the relative radii 




0.30174 


±0.00005 


0.30267 


±0.00039 


0.30295 


±0.00051 


Ratio of the radii 


k 


0.9291 


±0.0081 


0.9292 


± 0.0083 


0.9295 


±0.0083 


Surface brightness ratio 


.1 


0.9889 


± 0.0027 


0.9904 


± 0.0027 


0.9908 


±0.0026 


Orbital inclination (deg.) 


i 


76.862 


±0.012 


76.821 


±0.016 


76.800 


± 0.022 


Eccentricity term 


ecosw 


-0.000017 


±0.000006 


-0.000016 


± 0.000006 


-0.000015 


±0.000006 


Eccentricity term 


es'ma) 


0.00178 


±0.00026 


0.00198 


± 0.00027 


0.00201 


±0.00026 


Orbital period (days) 


Porb 


3.96004673 


±0.00000015 


3.96004673 


±0.00000015 


3.96004673 


±0.00000015 


Time of primary minimum 




52827.195695 


±0.000011 


52827.195694 


±0.000011 


52827.195693 


±0.000011 


Linear LD coefficient 


11 


0.380 


±0.017 


0.049 


±0.138 


0.191 


±0.062 


Non-linear LD coefficient 


V 






0.597 


±0.223 


0.273 


±0.098 


Relative radius of star A 


Th 


0.15642 


± 0.00065 


0.15689 


± 0.00066 


0.15701 


± 0.00068 


Relative radius of star B 


f B 


0.14533 


± 0.00075 


0.14578 


± 0.00073 


0.14594 


± 0.00076 


Light ratio 


tfi/tA 


0.853 


±0.016 


0.854 


±0.016 


0.855 


±0.016 


Orbital eccentricity 


e 


0.00178 


±0.00026 


0.00198 


±0.00027 


0.00201 


±0.00026 


Periastron longitude (deg.) 
Reduced x 2 


CO 


90.55 


±0.22 


90.46 


±0.20 


90.43 


±0.19 


y 1 


1.332 


1.328 


1.327 



rameters to avoid using theoretical coefficients. Because the two 
stars are very similar we also investigated solutions using the 
same LD coefficients for both stars. 



3.3. Light curve solution 

When obtaining the best fit to the available data we included 
the wire light curve (30015 datapoints), the Hipparcos H p light 
curve (131 datapoints), six photoelectrically-measured times of 
minimum light and the spectroscopic light ratio obtained by 
NJ94. The wire observations themselves have a root-mean- 
squared residual of 0.3 mmag per point. We found that the or- 
bital eccentricity is small but significant. When third light was 
included in the photometric solution we found that only a small 
and ill-defined amount was needed, resulting in very minor 
changes in other parameters. We have therefore assumed that 
the light curve contains no contribution in addition to those 
from the two components of the dEB. A mass ratio o f 0.97 was 
adopted from the spectroscopic study of Smith ( 1948), and grav- 
ity brightening exponents were set to 1.0 dClaretl fl998l) : large 
changes in these parameters have a negligible effect on the light 
curve as the stars are well separated in their orbit. 



When fitting the data we optimised + r%, k, the LD coef- 
ficients, the orbital inclination (/), the ratio of the central surface 
brightnesses of the two stars (/), the orbital period (P) and the 
reference time of minimum light (rMini). Orbital eccentricity, e, 
and the longitude of periastron, co, were included in the usual 
ebop way, by fitting for the quantities e cos to and e sin a> which 
have a more d irect dependence on the shape of t he light curve 
than e and a> ( lEtzell 1 1 975h iPopper & Etzellll98ll) . The out-of- 
eclipse magnitude of the system was also included as a fitted 
parameter. The wire data do not have a well-defined apparent 
magnitude zeropoint so this parameter has only a limited mean- 
ing and will not be discussed further. 

Our initial attempts to include the spectroscopic light ratio 
from NJ94 were only partially successful. However, this prob- 
lem disappeared completely if we included the reflection co- 
efficients as free parameters rather than calculating them an- 
alytically from the system geometry, ebop's simple analytical 
treatment of the reflectio n effect uses the un iformly-illuminated- 
hemisphere approach of Binnendiik ( 1974) and despite provid- 
ing an apparently good fit to the data it appears to be inadequate 
for observations of the quality of the wire light curve. This prob- 
lem can be avoided by including the reflection coefficients as pa- 
rameters of the fit rather than calculating them from the geome- 
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Table 5. Comparison between the measured LD coefficients for (3 Aur and predictions calculated using theoretical model atmo- 
spheres. Bilinear interpolation was used to interpolate the tabulated theoretical coefficients to the effective temperatures and surface 
gravities for ySAur given by NJ94. A solar metal abundance and a microturbulence velocity of 2kms~' have been adopted and 
results are quoted for the Jo hnson V and Cousins R bandpasses. 

Refere nces: V an Hamme is Van Hamme ( 1993 ); Diaz-Cordoyes is iDfaz-Cordoves et all {[995) and IClaret et alj (|T995); Claret is 



IClaretl d2000al) and Claret & Hauschildt is lClaret & Hauschildtl 12003). 



LD coefficient 


Star 


/JA 


ur 


Van Hamme 


Dfaz-Cordoves 


Claret (atlas9) 


Claret (phoenix) 


Claret & Hauschildt 










Vj 


Rc 


Vj 


Rc 


Vj 


Rc 


Vj 


Rc 


Vj 


Rc 


Linear law 


A 


0.382 ± 


0.006 


0.437 


0.367 


0.525 


0.418 


0.536 


0.542 


0.534 


0.454 


0.489 


0.406 


Linear law 


B 


0.379 ± 


0.006 


0.445 


0.372 


0.532 


0.402 


0.451 


0.455 


0.541 


0.459 


0.488 


0.404 


Sqrt (linear) 


A 


-0.06 ± 


0.15 


0.005 


-0.081 


0.016 


-0.019 


0.035 


-0.005 


0.008 


-0.044 


-0.245 


-0.240 


Sqrt (linear) 


B 


-0.18± 


0.15 


0.015 


-0.003 


0.028 


-0.017 


0.046 


0.001 


0.017 


-0.041 


-0.186 


-0.183 


Sqrt (nonlinear) 


A 


0.70 + 


0.23 


0.720 


0.625 


0.722 


0.598 


0.664 


0.603 


0.695 


0.658 


1.128 


0.993 


Sqrt (nonlinear) 


B 


0.88 ± 


0.24 


0.717 


0.625 


0.714 


0.599 


0.657 


0.602 


0.692 


0.661 


1.036 


0.902 


Quad (linear) 


A 


0.21 ± 


0.06 






0.273 


0.204 


0.208 


0.159 


0.203 


0.146 


0.242 


0.192 


Quad (linear) 


B 


0.19± 


0.06 






0.282 


0.196 


0.217 


0.164 


0.211 


0.150 


0.258 


0.206 


Quad (nonlin.) 


A 


0.25 + 


0.10 






0.328 


0.279 


0.387 


0.345 


0.401 


0.372 


0.366 


0.318 


Quad (nonlin.) 


B 


0.28 ± 


0.09 






0.325 


0.268 


0.383 


0.344 


0.400 


0.374 


0.340 


0.294 



try of the system. Their values are very well determined as they 
depend on the amplitude of the light variation outside eclipse, 
and this approach allows us to obtain an excellent fit to the light 
curve and spectroscopic light ratio. 

Measurement errors for the wire data are not available, so 
once this solution was obtained we calculated the appropriate 
measurement errors based on the residuals of the fit. We then 
binned these data by a factor of six, reducing the number of dat- 
apoints to 5027. In Paper I we noticed that a similar procedure 
during the analysis of if/ Cen resulted in fits which had a reduced 
chi-squared (x 2 [ed ) of 1.15 rather than 1.0, meaning that there 
was a non-random contribution to the residuals of the fit. We at- 
tributed this to a small systematic contribution to the photomet- 
ric errors in the wire data and multiplied the final uncertainties 

on each parameter by -^ 2 red to take this into account. Here the 

same effect is present, and the x 2 m & values are in the region of 
1.33. Minor systematic residuals in the fit might also contribute 
to^f 2 red , and in this case are possibly attributable to the approxi- 
mation used for the reflection effect. It must be remembered that 
the size of this effect is extremely small and would be essentially 
undetectable in a dataset with photometric errors greater than a 
millimagnitude. 

Photometric solutions for the assembled data (using the 
binned wire light curve) and different LD laws are given in 
Table|4] both with and without forcing the coefficients for the 
two stars to be equal. Solutions with linear LD have x 2 kA ~ 1-332 
whereas those using non-linear LD have x 2 rei ~ 1.328. The 
change in x 1 is about 20<x. According to an F-ratio test this dif- 
ference is significant beyond the 99.9% level, and we prefer so- 
lutions with non-linear LD (which have a lower ^ 2 red ) as we a 
priori expect these to be the most physically realistic. 

It is pleasing to discover that comparable solutions with the 
two different non-linear LD laws are overall in very good agree- 
ment. The greatest variation between all the solutions is seen 
for ta + t"b, which is clearly larger when non-linear LD is used 
compared to linear LD. This is qualitatively as ex pected (see 
Sectio n |3.2.1| >. We confirm the finding of Paper I and lLacv et al.1 
(12001) that using non-linear LD for dEB s causes the deri ved 
fractional stellar radii to be larger. Whilst IL~acv et al.l (120051) and 
Paper I found that the increase was or the order of 0.2%, here we 
find it is 0.3-0.5%. This can easily be attributed to the shallow- 
ness of the eclipses shown by (3 Aur. We also note that including 
non-linear LD coefficients causes the best-fitting orbital inclina- 



tion to decrease slightly, by 0.06°. This is substantially larger 
than the random errors found using Monte Carlo simulations. 

A comparison of the measured LD coefficients for j3 Aur with 
those calculated using theoretical stellar atmosphere models is 
given in Table|5] We have adopted the effective temperatures and 
surface gravities for/? Aur from NJ94. It is possible that the val- 
ues we obtain for the fitted LD coefficients are slightly affected 
by the ebop approximations to the stellar shapes and the reflec- 
tion effect, but this bias should be negligible. The interpretation 
of Table|5]is hindered because the photometric passband of the 
wire star tracker is not precisely known, but the theoretical pre- 
dictions are generally in line with the observed coefficients. The 
linear LD coefficients are a little smaller than most of the pre- 
dicte d ones, an effect which w as also seen in Paper I for t/f Cen 
and lSouthworth etaTJ ([2004b) for V453 Cyg. Further investiga- 
tion would be useful to see if this is an emerging pattern, and this 
effect should be borne in mind when modelling the light curves 
of dEBs using the linear LD law. 

3.4. Error analysis 

There are two separate contributions to the uncertainties of 
the derived photometric parameters: statistical and systematic. 
Statistical errors result from the measurement errors in the data 
and can be assessed using Monte Carlo simulations. Systematic 
errors may be caused by approximations in the model used to fit 
the light curve and can be estimated by using a different model 
to fit the data independently. Here we will do both. 

3.4.1. Statistical errors 

The formal errors which are often derived from the covari- 
ance matrices of the best fit to the light curve of a dEB can 
signi ficantly underestimate the true uncertainties in the solu- 
tion dPopperl Il984t iPopper et all Il986t [M aceroni & Rucinski 
Il997b ISouthworfh et alj l2005at ISouthwort h & Claused 120071 
This is because they do not fully take into account the 
correlations between differ ent photometric parameters (see 
IWilson & Van Hamrriel [2004). A good way of determining ro- 
bust statistical unc ertainties is by using Monte Carlo simulations 
(Pres s et al.l 19921 p. 684). T his has been implemented in jktebop 
(Southworth et al. 2004b c) and found to give results which are 
in excellent agreement with independent parameters determined 
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0.3015 

0.90 0.92 0.94 0.96 




0.90 0.92 0.94 0.96 0.9 



0.154 0.156 0.158 



-0.6-0.4-0.2 0.0 0.2 0.4 0.6 



-1.0-0.8-0.6-0.4-0.2 0.0 0.2 0.4 





0.37 0.38 0.39 0.40 




-0.6-0.4-0.2 0.0 0.2 0.4 0.6 



0.0 0.5 1.0 1.5 



0.149 
0.148 
0.147 



Quod- 







0.6 [ 
0.5 - 




0.154 0.156 0.158 



0.0 0.1 0.2 0.3 0.4 0.5 



0.0 0.1 0.2 0.3 0.4 



0.0 0.1 0.2 0.3 0.4 0.5 -0.10.0 0.1 0.2 0.3 0.4 0.5 0.6 



Fig. 3. Spreads in the values of different photometric parameters from the Monte Carlo simulations. 5000 simulations were done for 
each set of results shown here. The top row contains results for linear LD, the middle row for square-root LD, and the bottom row 
for quadratic LD. 



from different light curves of the same dEB (ISouthworth et al.l 
I2004bll2005cb . 

The uncertainties for the photometric parameters given in 
this work have been found using Monte Carlo simulations (5000 
per fit). We have quoted 1 cr (68.3%) confidence intervals for ev- 
ery quantity which has a value derived from each Monte Carlo 
simulation, except ^ 2 red . To reduce the amount of computing 
time required we used wire light curve after binning to reduce 
the number of points by a factor of six. 

The results of the Monte Carlo simulations are also excellent 
for searching for correlations between parameters. In this case 
the light curve is generally very well defined except for the ratio 
of the radii, k, as previously discussed. In Fig.[3]are shown plots 
of the spreads of values of various parameters, concentrating on 
those which are important or otherwise interesting. The poorly 
defined k means that r& and r% are quite correlated - as ta + r% 
is well defined their values are very strongly dependent on the 
spectroscopic light ratio. 

The LD coefficients are well defined in the case of the linear 
law, but it is clear that for the non-linear laws the two coefficients 
for each star are hugely correlated, which is not surprising. This 
seems to be worse for the square-root law than for the quadratic 
law, although this does not cause the overall light curve solu- 
tions for the two laws to be different either in absolute value or 
uncertainty. A much weaker correlation is seen between com- 
parable coefficients for the two stars. This shows that in gen- 
eral non-linear LD coefficients can only be measured using light 



curves of a comparable quality to the wire observations of (3 Aur. 
However, we would expect that dEBs with much deeper or total 
eclipses would give comparatively better results. A small clus- 
tering of points for the Monte Carlo simulations using quadratic 
LD is visible at Va ~ 0.0. This is due to a slight numerical prob- 
lem when calculating symmetrical partial derivatives. It is suffi- 
ciently far from the best-fitting value to have no effect except to 
slightly increase the derived uncertainties. 

3.4.2. Systematic errors 

ebop is suitable only for the analysis of well-detached eclips- 
ing binaries because it uses quite a simple approximation to the 
shapes of deformed stars. This approximation may be checked 
by fitting a light curve with a more complex code which 
uses Roche geometry to approximate the surfaces of the stars. 
We have theref ore studied the wire dataset using the Wilson- 
Devinney code dWilson & Devinnevll 19711 119731: Iwilsonll 19791) 
in its 2003 version (hereafter wd2003). 

As wd2003 does not have the option of directly fitting times 
of minimum light, we chose to consider only the (binned) wire 
data with the orbital ephemeris fixed to that obtained during our 
jktebop analysis. The square-root LD law was use d and the co- 
efficie nts were interpolated from the tabulations bv lVan Ha mme 
(1 19931) assuming a Cousins R passband. The linear coefficients 
were included as fitted parameters but the non-linear ones were 
not. The effective temperatures were fixed at the values given by 
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Table 6. Parameters of our solution of the wire light curve using 
wd2003. The optimised parameters are given in bold type and 
the fixed parameters in normal type. The rotation rates of the 
stars are given i n units of the synchronous v elocity. 
References: m lSmithl (11948b : (2) NJ94; (3) IClaret](ll998l) . 



Parameter Star A Star B Ref. 



1 



2 



3 



r (volume) 


0.15717 


0.14550 


r (point) 


0.15648 


0.14495 


r (pole) 


0.15807 


0.14621 


r (side) 


0.15707 


0.14541 


r (back) 


0.15787 


0.14607 



NJ94 and the wd2003 mode was used so the light contribu- 
tions of the two stars were not connected to these quantities. To 
force the solution to reproduce the spectroscopic light ratio from 
NJ94 we fixed the light contribution from star A to an appropri- 
ate value. This rather brute-force approach was needed because 
wd2003 does not have a facility to fix the light ratio directly. 
Due to this approach, the formal errors of the fitted parameters 
(which come from the covariance matrix) are unrealistic so we 
will not quote them. 

During the initial so lutions (with an integration a ccuracy of 
Nl = N2 = 30 points: IWilson & Van Hammdl2004 it was not 
possible to get a good fit to the light curve. The problem arose 
because the amplitude of the light variation outside eclipses was 
being underestimated, resulting in the profile of the eclipse be- 
ing very poorly matched. We were able to solve the problem 
by using the maximum integration accuracy (Nl = N2 = 60), 
including the bolometric albedos as fitting parameters, and ex- 
ecuting a large number of iterations. Due to the high integra- 
tion accuracy and the eccentric orbits required to fit (3 Aur, using 
the detailed reflection effect in wd2003 would have been pro- 
hibitively expensive in terms of calculation time. We therefore 
used the simple reflection effect, but preliminary analyses indi- 
cated that this makes a negligible difference to our results. Our 
final wd2003 solution has slightly larger residuals than for the 
jktebop solution, but the preliminary analyses have shown that 
this does not have a significant effect on the final parameters. 
Because of this our wd2003 analysis should not be regarded as 
definitive, and we hope to improve it by implementing a more 
modern minimisation algorithm into wd2003 (see discussion in 
Paper I). 

The final wd2003 solution is given in Table[6] and is in ex- 
tremely good agreement with the jktebop results. The great- 
est difference between the jktebop and wd2003 solutions is for 
the orbital inclination, which is 0.1° (4.6<x) hi gher in wd2003 
than ebop. A similar result has been found by Pop per & Etzell 



Table 7. Final photometric parameters for/? Aur from fitting the 
wire and Hipparcos light curves, the spectroscopic light ratio of 
NJ94 and several times of minimum light. The quadratic LD law 
was used with independently-optimised LD coefficients for the 
two stars. The upper part of the table contains directly optimised 
parameters and the lower part contains dependent parameters. 

The uncertainties have been multiplied by Jx 2 le ,\- 



Parameter 


Symbol 


Value 


Sum of relative radii 


r A + r B 


0.30290 


± 0.00054 


Ratio of radii 


k 


0.9300 


± 0.0097 


Orbital inclination (deg.) 


i 


76.80 


±0.10 


Eccentricity term 


ecosai 


-0.00002 


±0.00001 


Eccentricity term 


es'ma) 


0.00183 


± 0.00040 


Linear LD coefficient 


"A 


0.210 


± 0.074 


Linear LD coefficient 


»B 


0.187 


± 0.069 


Non-linear LD coefficient 


VA 


0.248 


±0.111 


Non-linear LD coefficient 


vb 


0.277 


±0.108 


Relative radius of star A 


f A 


0.15694 


±0.00081 


Relative radius of star B 


rB 


0.14595 


± 0.00082 


Light ratio 




0.854 


±0.016 


Orbital eccentricity 


e 


0.00183 


± 0.00040 


Longitude of periastron (deg) 


CO 


90.47 


±0.26 



d 19811) and lAndersen et al.l i 19931) when comparing ebop and 
wink dWoodll 197 111 1973b . and attributed to the treatment of the 
stellar shapes in ebop. However, the size of this effect is not 
enough to cause significant problems. This systematic error has 
been included in further analyses. For the other parameters there 
is no reason to believe that systematic errors are significant so 
we have retained their statistical errors as the final uncertainties. 

In the wd2003 analysis the bolometric albedos (Tabled con- 
verged to values near the canonical 0.5 for convective atmo- 
sph eres rather than the expected 1 .0 fo r radiative atmospheres 
(see lRafert & TwiggHl98CtlClaretll200lb . This is surprising be- 
cause the effective temperatures of the components of /?Aur, 
9350 K a nd 9200 K (NJ94) are appropriate for radiative atmo- 
spheres dClaretl 120005 iHurlev et all 120021) . We are also not 
aware of any similar cases in the literature. However, as we 
do not have realistic uncertainties for the best-fitting bolomet- 
ric albdeo values (see above), it is not clear how significant this 
discrepancy is. Further investigations of this effect will be un- 
dertaken in the future. 

3.5. Final light curve parameters 

It is very pleasing to discover that the six solutions listed in 
Table|4] give very similar results. The two with non-linear LD 
and separate coefficients for the two stars have the lowest ;f 2 red - 
Of these two the quadratic LD law solution has a marginally bet- 
ter X 2 K d an< ^ gi yes tne best-defined LD coefficients, so we will 
adopt this as our final solution. The two solutions differ by less 
than 0.1% in r& and r^. The quadratic-LD fit is shown in Fig. [2] 
and the times of minimum light are given in Table|3] We find that 
the ephemeris of the system is 

r Mi „i(HJD) = 2452 827.195693(12) + 3.96004673(17)£ 

where E is the cycle number and the quantities in parentheses 
indicate the uncertainty in the final digit of the preceding num- 
ber. The quoted uncertainty for the constant term corresponds 
to a probably optimistic 1.0 s, which is of a similar magnitude 
to small systematic effects such as the difference between helio- 
centric and barycentric Julian Date. 



Integration accuracy (points) 

Time of primary min. (HJD) 

Orbital period (day) 

Phase shift 

Mass ratio 

Orbital inclination 

Orbital eccentricity 

Longitude of periastron (deg) 

Mean efective temperatures (K) 

Potentials 

Rotation rates 

Bolometric albedo 

Gravity brightening exp. 



Light output 
Third light 



60 60 
2453827.19569 
3.9600467300 
-0.000001 
0.964 
76.909 
0.00203 
90.30 
9350 9200 
7.3309 7.5704 



1.0 
0.590 

1.0 
0.109 

0.624 
6.775 



1.0 
0.558 

1.0 
0.102 
0.624 
5.801 



0.0 
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Table 8. The physical parameters of the component stars of 
B Aur derived fro m the jktebop analysis and the circular spec- 
troscopic orbit of Smith (1 19481) . The effective temperatures are 
from NJ94. 

* Calculated ass uming L = 3.826xl0 26 W and M bo iO = 4.75 
dZombeckll990l) . 



Parameter 



Star A 



StarB 



Orbital separation ( R Q ) 




17.60 


±0.06 




Mass (M G ) 


2.376 


± 0.027 


2.291 


± 0.027 


Radius (R ) 


2.762 


+ 0.017 


2.568 


±0.017 


log(g/cm s~ 2 ) 


3.932 


±0.005 


3.979 


±0.005 


V sync h (kms -1 ) 


35.3 


±0.2 


32.8 


±0.2 


log Teg (K) 


3.971 


±0.009 


3.964 


± 0.009 


log(L/L Q ) * 


1.721 


±0.038 


1.630 


±0.038 


M bol * 


0.448 


± 0.094 


0.676 


±0.095 



The uncertainties in Table|7]come fr om th e Monte Carlo sim- 
ulations and have been multiplied by Jx*~ d = 1.152 to account 

for the small systematic errors detected in Section l3~3l A minor 
systematic error was found in the orbital inclination, and its un- 
certainty in TableQhas been increased to reflect this. 

The final photometric parameters are m ore accurate than and 
in good agreement with those of NJ94. iHummel et all d 1995b 
found an orbital inclination of i = 76.0° ± 0.4° from interfero- 
metric observations - this is 0.8° (2.0<x) lower than our result. 

4. Absolute dimensions and comparison with 
theoretical models 

The fundamental physical properties of the components of B Aur 
(TableO have been calculated from the final results of the jkte- 
bop analys is (Table[7]l, the circular spectroscopic orbit given by 
Smith (1948U and the effective temperatures given by NJ94. 
The masses are known to an accuracy of 1 .2% and the radii to 
0.5%, and the stars are slightly evolved with surface gravities of 
just below 4.0 (cgs units). 

The distance to BAm can be calculated in several ways 
(Table|9]l. The most direct of the se are the trigonometri c parallax 
measurements from Hipparcos dPerrvman et al.|[T997l) . and the 
orbital parallax obtained from the linear separation of the star s 
(Table[S]l and the orbital angular size from Hum mel et al.l (fl 995 ) . 
These two measurements are in excellent agreement. 

The distance to BAm can also be calculated us- 
ing surface brightness c alibrations or bolometric corrections 
dSoufhworth et al.l l2005a). both of which require reliable appar- 
ent magnitudes. We have corrected the mean Hipparcos Bj and 
Vt magnitudes to remove the effect of points wit hin eclipse, and 
supplemented these with JK magnitudes from iJohnsonl ([1964) 
adopting estimated uncertain ties of 0.02 mag for the J K data. 
Applying the calibration of Moo n & Dworetskvl d 1985b to the 
Stromgren indices given by NJ94, we find that the reddening is 
negligible. This is unsurprising given the closeness of B Aur to 
the Sun (25 pc) so we adopt E B - V = 0.00 ± 0.01. 

We have estimated the distance to BAm by th e surfac e 
brightness methocQ presented by [Southworth et al.l d2005a). 



3 ISmithl ( 119481) gives probable errors for the orbital elements he de- 
rives. Thes e must be mu ltiplied by a factor of 1.48 to give standard 
errors (e.g. iGriffinl 120051) , so the velocity amplitudes used here are 
K A = 107.46 ±0.58 km s~' and Kg = 1 1 1.49 ± 0.55 km s -1 . 

4 fortran77 code for calculating the absolute dimensions 
and various distance measurements for a dEB is available at 
http : //www . astro . keele . ac . uk/~ j kt/ 



which uses calibr ations of surface br ightness against effec- 
tive temperature dKervella et al.l 120041) . We have also cal- 
culat ed the distance using th eoretica l bolometric correc tions 
(see ISouthworth et all l2005bl) from iBessell et al] d!998l) and 
Girar di et al.l d2002l) a nd empirical b olometric corrections from 
Code et all d 19761) and [Flower! d 19961) . The good news is that the 
agreement between all four different methods of distance deter- 
mination is excellent. Whilst the two parallax distances are the 
most reliable because they are empirical direct measurements, 
formally the most accurate results are obtained using bolometric 
corrections or surface brightness relations in the K band, where 
uncertainties in effective temperature and interstellar extinction 
have the smallest effect. 

We have compared the predictions of several sets of 
theoretical stellar evolutionary models to the proper t ies of 
BAm. These sets a re the Granada models dClaredl 1995L Il997h 
[ Claret & Gimenez 1995), the Pad ova model s ( Girardi et al.l 
2000), and the Cambridge models dPols et alJll998l) . In each 
case a good fit can be found to the masses, radii and effective 
temperatures of jSAur for an approximately solar metal abun- 
dance (Z = 0.02) and an age of between 450 and 500 Myr. For 
a change in Z of 0.01 the quality of the fit is slightly decreased 
and the age changes by about 50 Myr, with a higher Z requiring 
a lower age. The properties of the two stars are too similar for 
a detailed comparison with theoretical predictions to be useful 
given the accuracy of the absolute dimensions we have estab- 
lished. An improved spectroscopic orbit will be needed before 
further conclusions can be made. 



5. Discussion and conclusions 

In this work we have presented a high-quality light curve of the 
detached eclipsing binary system jSAurigae, containing 30015 
datapoints with a point-to-point scatter of less than 0.3 mmag, 
obtained using the star tracker on board the wire satellite. Light 
curves of this quality will become available in the future for 
dEBs (from space missions such as Kepler and CoRoT), and it is 
important to ensure that the accuracy of light curve models will 
not simply be overwhelmed by these datasets. 

The data have been analysed using a version of the ebop code 
which has been extended to include non-linear limb darkening 
laws and the direct inclusion of observed times of minimum light 
and spectroscopic light ratios in the light curve model fit. Robust 
uncertainties have been obtained using extensive Monte Carlo 
simulations. The eclipses exhibited by B Aur are partial and quite 
shallow, so the ratio of the stellar radii is poorly defined by the 
photometric data. Including the accurate spectroscopic light ratio 
obtained by NJ94 causes the uncertainties in the stellar radii to 
drop by a remarkable amount, from 5% to 0.5%. 

The importance of including non-linear limb darkening laws 
in the analyses of dEB light curve has been investigated using the 
quadratic and square-root laws. An F-ratio test indicates that the 
resulting improvement in the fit is highly significant, and we find 
that the measured radii of the stars are 0.4% greater when using 
non-linear LD compared to the simple linear law, in agreement 
with previous studies. This bias is large enough to be important 
and we recommend that non-linear LD laws are used in photo- 
metric analyses of the light curves of dEBs. B Aur is not the best 
system to test this because its eclipses are shallow, so it is possi- 
ble that similar conclusions could be reached using lower-quality 
data of a deeply or totally eclipsing system. 

The LD coefficients obtained for B Aur for the linear, square- 
root and quadratic laws are in reasonable agreement with theo- 
retical predictions from model atmosphere calculations. We find 
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Method 




Distance (pc) 






Hipparcos parallax 




25.2 


±0.5 






Orbital parallax 




24.8 


±0.8 






B 


V 


J 


K 




Surface brightness calibrations 25.5 ± 1.8 


24.8 


±1.6 


24.8 ±0.5 


25.0 ± 


0.4 


Bolometric corrections (Bessel et al.) 


24.7 


±0.6 




24.8 ± 


0.3 


Bolometric corrections (Girardi et al.) 25.5 ± 0.8 


24.7 


±0.6 


24.5 ± 0.4 


24.8 ± 


0.3 


Empirical bolometric corrections (Code et al.) 


24.4 


±0.8 








Empirical bolometric corrections (Flower) 


24.9 


±0.9 









that the coefficients of bi-parametric laws are strongly correlated 
with each other, and that this will prevent their simultaneous 
measurement except for datasets of comparable quality to the 
wire light curve. The quadratic LD coefficients are better deter- 
mined than the square-root coefficients, which suggests that the 
latter law may actully be preferable as changes in the coefficients 
have less effect on the other parameters of the fit. The light curve 
fits for all three LD laws are in good agreement, so the choice of 
non-linear law is not important. 

In order to detect any systematic errors in the jktebop fit to 
the data, we have also fitted the wire light curve using the 2003 
version of the Wilson-Devinney code. We were able to obtain 
a reasonable fit but only by using the highest available inte- 
gration accuracy, bolometric albedos appropriate to convective 
rather than radiative envelopes, and a large number of iterations. 
The results agree well with the ebop solution, expect for a dif- 
ference of 0.1° in orbital inclination, which indicates that any 
systematic errors are too small to be important. 

The physical properties of the component stars of f3 Aur have 
been calculated from the results of the jktebop analysis and the 
spectroscopic orbit of Smith ( 1948), resulting in radii known to 
an accuracy of about 0.5% and masses to an accuracy of 1.2%. 
A good agreement between these properties and the predictions 
of theoretical stellar evolutionary models is found for an ap- 
proximately solar chemical composition and an age of 450 to 
500 Myr. This comparison is strongly dependent on the derived 
ratio of the stellar radii, which is itself determined primarily by 
the spectroscopic light ratio from NJ94. We will be obtaining 
new CCD spectroscopic observations of /3 Aur in order to calcu- 
late a more accurate spectroscopic orbit and light ratio, with the 
aim of measuring the masses as well as the radii of the compo- 
nents of /J Aur to accuracies of 0.5%. 

This analysis has shown that jktebop can provide an excel- 
lent fit to the high-precision wire light curve of the bright eclips- 
ing system fi Aur. Solutions using wd2003 indicated that there 
were no significant systematic effects in the jktebop solution. 
wd2003 was able to give a reasonable fit the the wire data, but at 
the expense of a large amount of computation time. We should 
therefore be able to derive highly reliable photometric solutions 
for the expected deluge of space-based eclipsing binary light 
curves from the Kepler and CoRoT missions. 
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